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NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 


TECHNICAL NOTE D-1870 

AN INVESTIGATION OF RESONANT, NONLINEAR, NONPLANAR 

FREE SURFACE OSCILLATIONS OF A FLUID* 
by R.E. Hutton** 

SUMMARY 

A theoretical investigation was conducted of the motion of fluid in 
a tank subjected to lateral harmonic vibration at a frequency in the 
neighborhood of the lowest resonant frequency of the mass of fluid. The 
investigation indicates that nonplanar fluid motion is due to a nonlinear 
coupling between fluid motions parallel and perpendicular to the plane 
of excitation, and that this coupling takes place through the free surface 
waves. These theoretical conclusions were experimentally confirmed 
using a cylindrical tank about 12 inches in diameter and partially filled 
with water to a depth of about 9 inches. The tank was translated along 
a diameter at selected frequencies in the immediate neighborhood of 
the natural frequency of small, free- surface oscillations. Three types 
of fluid motion were observed: stable planar, stable nonplanar (rotary), 

and unstable (swirling). In stable planar motion, the fluid moved har- 
monically with peak wave height and one stationary nodal diameter per- 
pendicular to the direction of excitation. In stable nonplanar motion, 
the fluid moved harmonically with a constant peak wave height and one 
nodal diameter that rotated at constant speed around the tank. In unstable 


Motion picture supplement L-77O has been prepared and is available on 
loan. A request card and a description of the film are included at the 
back of this document. 

** 

The author wishes to express his appreciation to Professor J. Miles, 
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motion the fluid never attained a steady-state harmonic response; peak 
wave height and nodal diameter rotation rate and direction changed 
constantly. 

The frequency regimes of the different types of motion which are 
possible in the neighborhood of the lowest natural frequency, f^, of 
small free- surface waves were as follows: 

1) Stable planar motion is possible except in a narrow 
frequency band roughly centered about f^. 

2) Stable nonplanar motion is possible in a band 
bounded below by f^. 

3) No stable motion, either planar or nonplanar, is 
possible in a finite frequency region just below f ^ . 


INTRODUCTION 

It has been widely observed that when a container of fluid having 
a free surface is subjected to a transverse harmonic vibration over some 
frequency ranges, the fluid does not necessarily respond with steady- 
state harmonic motion in which there is one stationary nodal diameter 
perpendicular to the direction of excitation. Rather, a wave may be set 
up which rotates around the tank harmonically or nonharmonically, even 
though the forcing motion is harmonic. The rotating wave phenomenon 
persists over a range of forcing frequencies centered about the natural 
frequency of small amplitude free surface waves. 

Specifically, the behaviour of a free- surface fluid in a container 
undergoing transverse harmonic vibration of increasing frequency is 
as follows. When the container is excited at a frequency appreciably 
below the lowest natural frequency, f^, of small, free-surface oscil- 
lations, the steady-state fluid motion is harmonic with a constant peak 
wave height and a single nodal diameter perpendicular to the direction 
of excitation. As the excitation frequency is increased, the wave height 
increases. At a frequency a little less than f^. 


the nodal diameter 
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begins to rotate at a nonsteady rate and with a varying peak wave height. 
This unstable swirling motion persists up to a frequency a little above 
the natural frequency, where once again the steady- state fluid motion 
is harmonic with constant peak wave height and a fixed nodal diameter 
perpendicular to the direction of excitation. An additional increase in 
excitation frequency from this stable point reduces the wave height 
until the cycle begins again as the next resonant frequency is 
approached.** 

In addition to harmonic planar motion and nonharmonic swirling 
motion, there is a third, which in this report is called harmonic non- 
planar motion. In the steady state, it is characterized by a constant 
peak wave height and a single nodal diameter which rotates at a constant 
rate. Harmonic nonplanar motion occurs in a frequency range bounded 
below by f ^ ^ . 

According to the linear approximations usually employed, a free 

surface fluid in a container undergoing transverse harmonic vibration 

should exhibit a steady- state planar harmonic motion at all frequencies 

except resonance. The fact that harmonic nonplanar and unstable 

2 

motions occur has lead investigators to suggest that nonlinear and 
viscous effects plan an important role. However, no extensive studies 
have been conducted to resolve this problem. The following work shows 
that the rotary and swirling motion can be accurately predicted in an 
inviscid liquid if the analysis includes the appropriate nonlinear effects. 

SYMBOLS 

tank radius 

coefficients in kinematic free- surface condition, 
Equation (A. 5) 

coefficients in kinematic free- surface condition, 
Equation (A. 6) 

Superscript numericals refer to references at the end of this report. 


a = 


a = 
n 


mn 
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b = coefficients in dynamic free surface condition, 
mn Equation (A. 8) 

c^ = amplitude perturbation coefficients, Equation (A. 42) 

d^^ = matrix elements defined in Equation (A. 44) 

f = first slosh mode generalized amplitude coefficients, 
n Equation (A. 23) 

^mn = natura ^ frequency of the mn'th sloshing mode, cps 
g = acceleration of gravity = 386 in/sec 
h = fluid depth 

i,j = unit vectors along x, y axes, Figure 1 
— * — ) — * 

i r> ig,i^ - unit vectors along r, 0, z axes, Figure 1 

kn» - parameters defined in Equation (A. 36) 

m, n = integer subscripts referring to the mn'th mode 

p = natural frequency of the mn'th sloshing mode, 
r mn , / 

rad / sec 

"q = fluid velocity vector 

r = fluid particle radial coordinate 

t = time 

x^ = tank displacement 

x^ = tank velocity 

A , B = velocity potential generalized coordinate 
mn mn 7 0 

A A A A 

A ,B ,C ,D = velocity potential generalized coordinates, 

mn mn nan mn. i-, . ■ /a ^ r\\ 

Equation (A. 30) 

B = free- surface boundary condition parameter. 


n 


Equation (A. 12) 


F^ = expansion coefficients defined by Equation (A. 38a) 
= functions defined by Equation (12) 
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H = function defined by Equation (13) 

A 

= parameters defined by Equation (A. 36) 

I = integrals defined by Equation (A. 36) 
n 

= integrals defined by Equations (A. 32) and (A. 36) 

J = Bessel function of the first kind of order m 
m 

K = parameter defined by Equation (A. 36) 

K q = parameters defined by Equation (A. 32) 

Kjq, K^q - parameters defined by Equation (A. 38b) 

Kj , = parameters defined by Equation (A. 38a) 

K^, = parameters defined by Equation (A. 41) 

= determinants defined by Equation (A. 47) and (A. 54) 

a mn - parameter defined in Equation (A. 36) 

y = velocity potential function steady- state amplitude, 
Equations (A. 39a) and (A. 40a) 

€ = tank velocity amplitude 

=. tank displacement amplitude 

t, = velocity potential function steady- state amplitude, 
Equation (A. 40a) 

t| = wave height 

9 = cylindrical coordinate system coordinate 

X = perturbation parameter in Equation (A. 42) 

, X = roots of perturbation characteristics equation 
/ 

\ - roots of J (\ a) = 0 

mn m mn 

v = transformed frequency, Equation (A. 19) 
t = transformed time, Equation (A. 19) 
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th 

k derivative of velocity potential function 
evaluation on z = 0, Equation (A. 4) 

angular forcing frequency of tank motion 

function defined in Equation (A. 1) 

total velocity potential function 

disturbance velocity potential function 

parameters defined by Equation (A. 31) 

disturbance potential functions defined by 
(A. 18), (A. 23) and (A. 27a, b) 



ANALYSIS 

The object of this investigation is to determine the response of a 
fluid contained in a cylindrical tank which is undergoing lateral harmonic 
vibration in the neighborhood of the lowest natural frequency of small, 
free surface waves. In this section the main steps in the analytical 
derivation are presented, with the algebraic details relegated to the 
appendix. 

Figure 1 illustrates the geometry and coordinate system which is 
used in the analysis. The tank is assumed to be forced in the fixed 
x o Z Q -plane only and this motion is denoted by x^(t). The r, 0, z coor- 
dinate system moves with the tank with the plane z = 0 coinciding with 
the quiescent free surface. 

Throughout the analysis it will be assumed that the fluid is 
incompressible, inviscid, initially irrotational and that the only body 


K = 

u> = 

r = 
$ = 
$ = 


Q n , n, = 

On 2n 

x , \\t = 

n n 
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force is that of gravity, which is oriented along the negative z-axis. 
Under these assumptions the fluid velocity vector c| can be written as 

q = V$ (1) 

where the velocity potential function $ is expressed in terms of the 
moving coordinates, r, 9, z. If $ is written as the sum of a disturbance 
potential function $ and a function accounting for the tank motion, i.e. 

$(r, 0, z,t) =x^ r cos 8 + $ (r, 0, z, t) , (2) 

then under the above assumptions the boundary value problem to be 

3 

considered can be stated as 


V 2 $(r, 0, z, t) 


$ =0 
z 


= o 



(3) 


$> = X, cos 0 

r 1 


on r = a 




and 


gr| + (v$ • V$) = - x i r cos 0 on z = t), 


(4) 


»>✓ 1 

T] + $ T) + — = $ 

’t r ‘r 2 0 1 0 z 

r 


on z 


(5) 


The dots above represent time derivatives and the subscripts 
t, r, 0, and z respectively. The first equation of (3) requires that the 
disturbance velocity potential function satisfy Laplace’s equation which 
is a linear partial differential equation. The dynamic condition, Equation 
(4), is Bernoulli’s equation with the pressure set equal to zero on the 
free surface z = q. The kinematic condition, Equation (5), is simply a 
statement that a fluid particle on the free surface has the same vertical 
velocity as the free surface. 
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It is to be noted that the only nonlinear character of this boundary 
value problem enters through the free- surface boundary conditions on 
z = rj . In a linear analysis, the nonlinear terms in Equation (4) and (5) 
are neglected under the assumption that the wave height and fluid velocities 
are small, and then, these linerized free- surface conditions are satisfied 
on the undisturbed free surface z = 0. However, when the tank is driven 
at or near the lowest resonant frequency f^, the wave height and fluid 
velocities are not small, and it is essential that the nonlinear terms in 
Equation (4) and (5) be taken into account. 

Since the free surface z = rj is an unknown in the problem, the 
free- surface boundary conditions are combined in the appendix into one 
boundary condition in which the wave height r\ has been eliminated. Thus, 
the boundary value problem to be considered can be restated as: 


f 0 rf r <a 

(r, 0, z, t) = 0 in < 0 I? 0 « 2tt 

z « n 

<J> = 0 on z - -h 

z 

% - 0 on r = a 

r 

and 




> ( 6 ) 


J 


B 1 + B 2 + B 3 + 0(t! 4 ) = 0 


(7) 


where the terms B . B , B , defined in the appendix, depend only upon 

1 2 3 ^ 

the velocity potential function and its derivaties, all evaluated on the 
undisturbed free- surface z = 0, and also upon the prescribed tank 
displacement x^(t). The tank displacement is taken as 

x.(t) = € sino.>t 

r o 
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with "small” and u> close to or equal to the lowest sloshing natural 
frequency p^. A steady state harmonic solution to boundary value 
problem (6) and (7) is posed in the form 


<$ = € 


^(r, t) cos u>t + x^(r, t) sin wtj 

^ [V r) + COi5 + sin 


+ € ijj^(r) cos 3cot + x^(r) sin 3cotj (8) 

where e = cot ^ is the peak tank velocity and the functions x^, ^ satisfy 
(6) identically with only *|/j and x^ depending upon time. The notation 
r means dependence upon r, 0, and z. A set of normal sloshing modes 
which satisfies (6) exactly is 


Ta (t) cos m9 + B (t) sin m0 
j mn mn 


J 


(X r) 
m mn 


cosh X (z+h) 
mn 

cosh \ h 
mn . 


where the J are Bessel functions of the first kind and of order 
m 

m = 0. 1, 2, • • • , and X are an infinite set of numbers for each m 

mn f 

obtained from the zeros of J (X a). The functions A (t) and B (t) 

m mn mn' mn' 

depend only upon time and are called the generalized coordinates of the 

mn'th mode. The natural frequency of small, free- surface oscillations 
in the mn'th mode is denoted by p . In this problem, in which the 
tank is harmonically translated at a frequency close to or at the lowest 
natural frequency associated with the mode, p^, g enera ^ zec l 

coordinates A^ and B^ dominate all other generalized coordinates. 

For this reason, the first order terms of (8) contain the first mode 
only and x^ and \\j^ are taken in the appendix as 
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vp = jf^(r) cos 0 + f 3 (r) s * n 


cosh 


J,(X,ir) 


^1 i( z+ h) 


1 11 cosh I X, , h 


rn h j 


x i = 


^(t) cos 0 + f^(t) sin 0: J^(\^r) 
where the transformation 


cosh X, . (z+h) 

■l_j__ J 

cosh X, . h 

L 11 1 


1 2/3 2 2 f. 2/3, 

T = -£ € Wt » p^=CO [ 1 - V € J 


(9) 


(10) 


has been introduced. The functions and x^ depend upon the modes 
corresponding to J q and only and their definitions are given in the 
appendix . 

In the appendix the trial solution (8) is introduced into the free- 

1/3 2/3 

surface conditon (7) and the terms of which « , € , and e are the 

1/3 

coefficients are set equal to zero. The term of which € is the 

coefficient depends upon the first J. mode only and vanishes identically. 

2/3 1 

The term of which e is the coefficient depends upon the generalized 
coordinate of the first mode and the infinite set of generalized coor- 
dinates of the J q and modes. With Fourier-Bessel techniques, this 
equation is satisfied by expressing the generalized coordinates of J q 
and modes in terms of the mode generalized coordinates. The 
term of which € is the coefficient depends upon the first Jj mode and 
coupling between the mode and J q and modes. 

The generalized coordinates of the J q and modes are eliminated 

from the last equation by using the relations (obtained by setting equal 

2/3 

to zero, the term of which e is the coefficient) between the generalized 
coordinates of the mode and the J q and modes. The resulting 
equation depends only upon the first Jj mode generalized coordinates. 
Next, the equation is satisfied in a Rayleigh-Ritz sense. That is, the 
equation is first multiplied by Jj cos 0 r dr d0 and integrated over the 
free surface and is then multiplied by sin 0 r dr d0 and again integrated 
over the free surface. 
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Both of these integrated equations have sin cot and cos ut terms. 
Requiring that the coefficients of each of the sin cot and cos cot terms 
vanish in both of these two integrated equations gives the nonlinear 
differential equations that the generalized coordinates fj, f^, iy f^» 
of the first mode must satisfy. This set of first order differential 
equations can be written in a compact form using Cartesian tensor 
summations and writing: 


df. 


It =0 i<V f 2' f 3- f 4> ' *= 


where 


C 1 =- H '2 


C 3 = -H. 4 


V H ’l 


°4 = H ' 3 


and 


H = F l f I + 7 vt ii + ? K l'W 2 -7 K 2< f 2 f 3 - f l f 4> 2 


( 11 ) 


( 12 ) 


(13) 


Subscripts following commas imply differentiation with respect to the 
corresponding f^ and repeated subscripts imply summation. The con- 
stants Fj, Kj , K which depend upon tank geometry only, are defined 
in the appendix. A •steady state harmonic solution of the boundary 
value problem expresed in Equations (6) and (7) corresponds to the 
zeros of the four equations indicated in (11). In the appendix it is 
found that there are two such solutions. The first is called planar 
motion and the solution is given by 


f l=Y i f 2 = f 3 = f 4 = 0 (14a) 

with 

^ -1 „ 2 
v = -Fj-y - Kjy 


(14b) 
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The second is called nonplanar motion and the solution is given by 

F 

f l=V , f 2 = f 3 = 0 , f 4 2 =b 2 = Y 2 + K l y" 1 (15a) 

2 

with 

v = -K 3 y _1 + K 4 y 2 • (15b) 

The constants K, and K are also defined in the appendix. In the 
j 4 

appendix these two solutions are examined to determine whether they 
are stable by imposing a slight perturbation from the steady state solution 
and investigating the subsequent motion. If the subsequent motion in- 
creases with time following the perturbation, the solution is called un- 
stable. The analytical investigation of stability is carried out in the 
appendix, wherein the perturbation 


f . ( t) = f^ + c. e^ T (16) 

is imposed. The f.^ denote the steady state solution given by (14) and 
(15) with the disturbance c. assumed small. Introducing (16) into the 
set (11) leads to an algebraic equation in X. Any root X of this equation 
with a positive real part indicates a motion that grows with time and the 
corresponding solution is called unstable. It is found that the solutions 
given by (14) and (15) are stable only over certain ranges of transformed 
forcing frequency v. Figure 2 summarizes the regions of stable and 
unstable planar motion. From the figure it is noted that the theory 
predicts that planar motion is not stable over a frequency band centered 
about v = 0 (or co = Pjj)- 


NUMERICAL EXAMPLE 

The following numerical calculations are based on the theoretical 
results for the parameters corresponding to the sloshing tests, which 


are 
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a = tank radius = 5. 938 inches , 

h = water depth = 8.907 inches . 

With the values \j^a= 1.84119 and Jj(kjja) = 0.581865, 
Equations (A. 24) and A. 38a) are evaluated to be 

Fj = 8.53992 , a n = 0.99205 

and 


Pjj = 10.897 rad/sec = 1.734 cps . 

In the evaluation of the coefficients and the infinite series terms 
in Equations (A. 35a, b) are approximated by their first five terms. To 
indicate what errors this finite series approximation might cause, the 
following facts are cited. In the calculations of Cr^ and 6^, which con- 
tribute to Kj and K^, the last three terms were only about 1 percent of 
the first two terms. 

The Bessel function parameters required to sum the finite series 
obtained from Reference 4 are presented in Table 1. 

The values in Table 1 are used to compute the quantities defined 
by Equations (A. 31), (A. 32), (A. 36), and (A. 38a, b) . The results are 
presented in Table 2. 

The amplitude- frequency relation for the planar motion is 


= -f^" 1 - KjV 2 = -8.5399Y -1 - 0.48528 x 10 _5 y 2 . (17) 


V 
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The lower limit boundaries between stable and unstable motions are 


computed from Equations (A. 49a, b) 


v ’ V*i 


= 95.82 




-0.1337 , 


The upper limit boundaries are computed from Equation (A. 50a, b): 


= -85.41 


< K z- K i»hc 


= 0.06459 


The amplitude-frequency relation for the nonplanar 


motion is 


v = -K 3 Y _1 +K 4 y 2 = -3.0235V" 1 + 4.0010 x 10" 6 y 2 (20) 


i Z = v 2 V" 1 = V 2 + 6.2305x 10 5 y _1 


Stability of the steady- state nonplanar solution is determined by 


examining the roots of Equation (A. 53), namely: 
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\ 4 + 


<m 3 + 


m 4 )\ 2 + 


m 5 m 6 


= 0 


( 22 ) 


where the coefficients are obtained from (A. 56) and are found to be 


+ M 
3 4 


2 / 3 . 3K 2 ‘ 2K ! F A 

2 — * 1 

4K. 


4K 2 S(y + 


7. 5148 x 10“ 10 v(y 3 + 3. 5699 x 10 5 ) , 


-1/3.113 


K 1 F 1 


M 5 M 6 = 4K, (K 2 - 2 Kl )F lY Iy +^iV ~ 7KJ 


= 2. 5677 x 10 " 9 y" 1 (Y 3 + 6 - 2305 xl 0 5 )(Y 3 + 3 - 7784 x 10 5 ). 

(23) 

Stability calculations for the nonplanar motion are summarized in 
Table 3. 


The computed results for this numerical example are presented 
in Figures 3 through 8. Figure 3 is a plot of Equation (17) showing the 
relation between y and v for the harmonic planar motion and indicating 
the corresponding ranges of stable and unstable motion. Figure 4 is 
a plot of Equation (20), showing the relation between y and v for the 
harmonic nonplanar motion; the ranges of stable and unstable motion 
are indicated, as well as the range in which the solution is imaginary. 
In Figure 5, the plot of t, versus v for this same nonplanar motion is 
derived from Equations (20) and (21). Figure 6 is a plot of y versus 
v, showing only the stable branches of Figures 3 and 4. Figure 7 shows 
the unstable regions for both the planar and nonplanar motions as 
functions of the peak tank velocity, €, and the tank displacement fre- 
quency, co/2tt. These curves are computed from the transformation 
equation 
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p n 2 = W 2 (l - v^/3, (24) 

for the value of p^ = 10.897 rad/ sec = 1.734 cps. An evaluation of 
Equation (24) requires the limiting values of v. The limits of v sepa- 
rating stable and unstable regions for the planar motion are shown in 
Figure 3 and are given in Equations (18) and (19) as v = -0. 1337 and 
v = 0.06459. The value of v separating stable and unstable regions 
for the only stable nonplanar motion, which is for y > 0, is shown in 
Figure 4 and is tabulated in Table 3 as v = -0.03027. Substituting 
these values into Equation (24) and solving for e as a function of w/2tt 
provides the data used to plot the curves in Figure 7. Figure 8 is a 
plot of the approximate planar motion stability boundary as a function 
of the dimensionless ratios f/f^ and € Q /a. These curves approximate 
the stability boundary for a tank in which the fluid depth is greater 
than the tank radius, in which case, the stability boundary is nearly 
independent of the fluid depth. The expressions used to construct this 
figure were developed in the following manner: First by referring to 

Equations (A. 36), (A.38a,b), (A. 49b) and (A. 50b) it is noted that the 

-1/3 

instability limits are approximately proportional to (ga) . That is 

f F l Y= /3 -1 -i( a '\ 2/3 , 

v u = ' 3K l(2lq-j “ 8 a I ■_ t ^_ 3 J = (gal 

”1 * < K 2 - K l > (^) V3 00 ( « a) ' I/3 

or, introducing the proportionality constants c^ and c g gives 

/ i-l/3 . . ,-1/3 

v u = c u (ga) and V i = C f (ga) 

For a = 5.938 and g = 386 the upper and lower stability limits are found 

to be v = -0. 1337 and v . = 0. 06459 so that c and c . become 
u i nl 



c = 0.852 
u 


c = -1.764 • 
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Using Equation (A. 19) the ratios of driving frequency to the lowest 
natural sloshing frequency can be written 



the upper and lower stability limits can be approximated by the 
expressions 



These were the equations used to develop Figure 8. 

Experimental Results and Comparison with Theoretical Predictions 

The sloshing tests for verification of the theoretical work of 
this investigation were conducted in the Space Technology Laboratories 
slosh test facility, shown in Figure 9. This facility consists of a 
rigid, lightweight tank platform supported so as to allow only a single 
translation degree of freedom. This is accomplished by long, flexure- 
ended positioning struts oriented to eliminate rotation about any axis 
and to make negligible any translation other than that along the lateral 
drive axis. This design results in a very rigid platform possessing a 
high resonant frequency compared with the slosh frequency being 
investigated. The drive mechanism and its platform are also very 
rigid. 
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The sinusoidal lateral oscillation of the slosh-tank platform is 
accomplished with a scotch-yoke mechanism connected to a heavy, 
variable- speed drive. The amount of eccentricity of the drive stud 
can be preset to produce a peak-to-peak stroke length up to 1 inch. 
While the drive motor is running, the drive stud can be moved rapidly 
from center (no stroke) to preset stroke position to start lateral exci- 
tation. Thus, the amount of excitation can be varied at any selected 
frequency. 

A cylindrical tank about 12 inches in diameter, partially filled 

with water to a depth of about 9 inches, was used for the tests. Testing 

was conducted in the following manner. The amplitude of tank motion, 

e , was set and was left constant for the entire test. Then the tank 
o 

was driven at a desired frequency and when steady- state conditions 
appeared to be established, the frequency and wave height were read 
and the fluid behavior noted. 

The driving frequency was measured by sensing the rotational 
speed of the scotch-yoke drive and recording this speed on a Berkeley 
freqeuncy meter. The tank peak-to-peak amplitude ( Ze q ) was meas- 
ured by placing a dial indicator beside the tank in the plane of tank 
motion and observing the extreme positions of the dial indicator needle. 
The wave height was measured by reading a 6- inch scale (graduated 
in hundredths of an inch) attached to the plexiglas tank in the plane of 
the tank motion and zeroed on the quiescent water level. The error 
of the three measurements is estimated to be as follows: drive fre- 
quency, a>/2Tr = ±0.004 ops; tank amplitude, € q = ±0.0005 inch; wave 
height was more difficult to read and the probable wave height error 
was ±0. 2 inch. ) 

Three types of fluid motion were noted: stable planar, stable 

nonplanar, and unstable. Stable planar motion is a steady- state fluid 
motion with a constant peak wave height and a stationary single nodal 
diameter perpendicular to the direction of excitation. Stable nonplanar 
motion is a steady- state rotary fluid motion with a constant peak wave 
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height and a single nodal diameter that rotates at a constant speed. 

This motion has the appearance of a surface wave traveling around 
the tank at a constant speed and in a single direction. Unstable motion 
is a swirling fluid motion that never attains a steady- state harmonic 
response; the peak wave height and nodal diameter rotation rate and 
direction continually change with time. The surface wave in unstable 
motion may build up and then decay; it may rotate first in one direction, 
then stop, and rotate the other way. At times, the wave may slosh 
for several cycles in a plane perpendicular to the plane in which the 
tank is being driven, then rotate around and slosh in the driven plane 
for several cycles. 

To locate the limits of the planar motion instability region, the 
first test reading was made at a frequency in the stable region well 
below the instability region. After the first reading was recorded 
the driving frequency was increased. When steady-state motions were 
again established, the wave height and fluid behavior were again re- 
corded. This process was continued until the frequency was raised to 
the point where planar motion became unstable, indicating the lower 
limit of the unstable planar region. With this boundary established, 
the driving frequency was increased above the natural frequency to a 
frequency in the upper stable planar region, and the steady- state wave 
height was recorded. This process was repeated, reducing the fre- 
quency each time, until the upper limit of the unstable region was 
located. 

After the planar motion exploration was completed, the stable 
nonplanar motion was investigated. This motion occurred primarily 
above the natural frequency and overlapped a frequency range in which 
stable planar motions had been observed. 

To verify further that both planar and nonplanar motions could be 
stable within a common frequency range, the following test was carried 
out. With the fluid at rest, the tank motion was started. When a 
steady- state condition' at a selected frequency was reached the fluid 



20 


was observed to assume stable planar motion. While the water main- 
tained this stable motion, a paddle was inserted into the tank and the 
motion was deliberately perturbed into a random behavior. When the 
paddle was removed, the water motion sometimes settled back into the 
original planar mode, and at other time became nonplanar . Observa- 
tions of the latter cases indicated that this nonplanar mode was indeed 
a stable motion which would continue until the driving motion was 
changed. This test thus demonstrated that both stable motions were 
possible at a single driving frequency. A possible conclusion from 
several repetitions of this experiment is that the nonplanar mode re- 
quires the addition of energy into the Bystem. 

Additional testing verified the existence of a finite region of un- 
stable motion below the natural frequency. In this region, no stable 
harmonic steady- state motions were found. 

The data taken during these tests are summarized in Figures 10, 
11, and 12. Figure 10 is a plot of faired curves through the observed 
planar motion peak wave height data for various tank driving fre- 
quencies. Figure 11 shows a faired curve through the test data of 
the nonplanar motion peak wave height for a driving amplitude of 0.032 
inches. Figure 12 is a plot of the data of Figure 10 when the frequency 
is transformed by the equation 



1/3 

and the wave height is divided by e . This scaling was suggested by 
the theoretical investigation. 


In Figure 13, the predicted stable and unstable planar and non- 
planar regions of Figure 7 are compared with the test data. Figure 13 
indicates that the predicted boundaries of the instability region agree 
closely with the experiment for small values of «. Poorer agreement 
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for the nonplanar motion stability boundary is not surprising because 
of the difficulty of visually detecting the point at which the fluid motion 
changes character; as the frequency is decreased, the transition is 
gradual from the steady- state harmonic nonplanar motion to the non- 
harmonic nonplanar motion existing in the instability region. 

CONCLUSIONS 

This investigation has attempted to provide an explanation for 
the motion of fluid in a partially filled tank when the tank is laterally 
vibrated over a frequency range centered at the lowest natural fre- 
quency of small free- surface oscillations. The theoretical studies 
outlined herein predicted the forcing frequency ranges on which there 
were stable steady- state harmonic planar and nonplanar fluid responses 
and the range of frequencies in which no harmonic motions are stable. 
These studies have shown that the mechanism that causes the swirling 
fluid motion in the unstable regime and the rotary harmonic motion 
in the nonplanar regime is a nonlinear coupling of fluid motions parallel 
with and perpendicular to the excitation plane. It is evident that this 
nonlinear coupling takes place through the free- surface waves and the 
observed fluid behavior can be predicted only if the higher-order terms 
in the free- surface dynamic and kinematic boundary conditions are 
included in the analysis. 

Space Technology Laboratories, Inc. 

Redondo Beach, California, November 1, 1962 
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APPENDIX 

1. Free Surface Boundary Condition 

Since the free- surface height r) is an unknown in the problem, 
it is desirable to replace the two free- surface conditions (4) and (5) by 
one equation that does not contain q. First, it will be convenient to 
eliminate the partial derivatives of q occurring in Equation (5) by use 
of Equation (4). Now, Equation (4) is of the form 

- gq = r[r, 0, q(r, 9, t), 0 = ? t + j fe 2 + \ 'Zq + S^+Xj r cos 0 , (A. 1) 

so that q is no longer an independent variable. Thus, the partial deriva- 
tives q , q , q to be obtained from Equation (A. 1) must reflect the 
t r 0 

fact that r also depends upon r, 0, and t through the free-surface param- 
eter q. From Equation (A. 1) it is determined that 


- gq t = r t + r n q t , 

- gq r = r r + r Ti q r . 

- gq 9 = r e + r Ti q 0 , 


■ (g + Tq)q t 

- (g +rq)q r 

- (g + Tq)q e 


r t = *tt + ®r*rt + “I®0®0t 

r 

r r = *rt + ®r*rr'“3*0 2 + 

r 

F 0 = ^0t + *r*r0 + 7Z®0*00 


rv fw i • • 

+ $ $ , + x.r cos 0 , 
z zt 1 

~T*0®0r + f z $ zr + *1 
r 

+ $ » - x, r sin 6 , 

z 0z 1 


N 

cos 0,> (A. 2) 


J 


I 
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where 


rv '“nj 1 n-w 

r = $ A + 9 9 + —,$-.$> + $ $ 

r) nt r r n 2 0 0 r| z t|,z 


Since T is equal to T for z - ri 
V z 


r n = F z = »zt + ®r*rz + 72*9»ez t »z® 


zz 


If Equation (5) is multiplied by -(g + F^) and Equations (A. 2) are used, 

Equation (5) reduces to 

_ «'*W ~ 2 ^ **** . ^ 2rNJ «%< 2^>» 

ff) 4* £<J> + Z<b & , 4 — =5$rk$rki. * ^$ $4.4'$ $ *$ $ 

w tt B z t rt 2 0 0t zt r rr z zz 

r 


+-?$ e 2 ?ee - -T^r^0 + '^r^0 $ r0 + ^Z*z%%9z + 2 *r*z^ 
r r r r 


rz 


- - Xj r cos 0 + 




$V. sin 0 - $ cos 0 
0 r 


)• 


on z= (A- 3) 

Since the potential functions must be evaluated on z = r\ t Equation 
(A. 3) depends upon q implicitly and Equation (4) depends upon t \ both 
implicitly and explicitly. Now the wave height r] can be eliminated 
between Equations (4) and (A. 3), if these two equations are first ex- 
panded in a Taylor series in r\ about z = 0. By introducing the notation 


5 $\(r, 0> z = 0, t) , 


(A. 4) 


where the subscript k represents any order of partial differentiation, 
the Taylor series expansions for Equations (A. 3) and (4) can be written 
in the forms 
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and 


2 3 

a Q + a^ + + a^ + . . . = 0 


b 0 + + ^2^ + ^ 3 ^^ + • * • = 0 » 


respectively, where 

a o = a oo + 

a l = a ll + 
a 2 = a 22 + 


• • • 


and where 

a oo = + + r *i cos 6 • 

a 01 = 2 ^r^rt + 72 Vet + 2 ^z^zt + ,3i l(^r COS 9 '7^9 8in e ) ' 

a 02 = ^rr + “T *0^00 + ^zz ' ~3^r^0 2 + 2 *rVrz 
r r 

+ “T * r * 0 *r 0 + “T 0 zV 0 z ’ 
r r 

1 


l ll 


= ***„ + g* 


ttz + g *zz ' a 22 " 7 ' ' ttzz 


+ g* zzz > * 


a 12 = 2 v rz^rt + Mrzt + “^0**01 +_ 7 *0*0zt + *zz*zt 

v r r 

+ *z*zzt) + ^](^rz 008 6 - 7 *»z sin 9 ) • 


00 = + ”1 cos 8 ’ b ll = *tz ’ b 22 = 2*t 


zz 


01 * z(»r 2 + jz *S + *z) ■ 


b 12 = *r*rz + “2 V 0 z + *z* 


ZZ 


r 

1 


9a mn , 1 9b mn 

a m+l,n+l (m + 1 )! 5z ’ m+l,n+l (m + 1 )! 5z 


(A. 5) 
(A. 6) 


+ a 02 * 

b 0 = 

o 

o 

+ b 01 ' 


+ a 13 , 

b i = * + b n 

+ b 12 ! 

V * (A. 7) 

+ a 24 * 

b 2 = 

b 22 

+ b 23 * , 



> (A. 8) 
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From Equations (A. 6) and (A. 7), it is evident that the 
potential functions are the 'same order as the wave height. This is 
readily seen by neglecting the products of q and <f> in Equation (A. 6) ; 
then, the first approximation becomes 


b 


00 




= 0 


or 

r\ t ~ — <l> t - (A. 9) 

2 3 

Hence, a term such as b^q in Equation (A. 6) is of order q . 

Solving Equation (A. 6) for q yields 


n 


B 7 -E 7 





0 + 0<T, 3 ) 


(A. 10) 


and upon substituting Equation (A. 10) into (A. 5), the result becomes 


a o + a i 




0(q 3 ) 


) + aZ \l ^ + ^ 


12 4 

+ 0(q*) = 0 


or 


a , b b n . 

l_0 +a 0 + 0 (ti 4 ) = 0 

1 b l 


(A. 11) 


Now, 


a I b 0 


<a l 1 + a 12 + a 13 ,<b 00 t hoi 1 


* + b n + b i3 


a ll b 00\ 


a l l b 01 + a 12 b 00 

V 6 


a l 1 b 00 b l 1^ 

2 J 

g / 


+ 0(q 4 ) 



26 


and 


a 2 b 0 _ ( a 22 * a 23 + a 24^ b 00 * b 01^ 
b l (g + b H + b j 2 ^ 


a 22 b 00 


1+ ou 4 ) 


t 


so that Equation (A. 11) becomes 

Bj + b 2 + b 3 + OU 4 ) = 0 


(A. 12) 


where 


B, = a 


00 * 


B, = a 


a ll b 00 


01 - g 


> (A. 13) 


B-, = a 


02 


a l 1 b 0 1 +a 12 b 00 , a ll b 00 b ll +a 22 b 00 
+ z 


g 


g 


In summary, the problem statement in terms of disturbance 
potential functions with the higher -order approximation for the free 
surface condition becomes 


r<a 

V 2 $(r, 9, z, t) = 0 in / 0^0 ^2 tt 

^-h<z<n , 


^ = 0 on z = -h , 

z 

= 0 on r = a , 
r 




) (A. 14) 


J 


B i + 


B 2 + 


1 3 + 0( T1 ) = 


on 


= 0 


(A. 15) 
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2. Steady-State Harmonic Solution 


In this section, a steady- state solution to boundary value problem 
^Equations (A. 14) and (A. 15)jis derived for the case where the tank 
displacement motion is 


x, = € sin cot 
1 o 


(A. 16) 


and the tank velocity is 


x i = 


€ COS COt 


€ = G0€ 


(A. 17) 


The solution, which is limited to the case where e is small and the 

driving frequency to is close to or equal to the first slosh mode fre- 

* 

quency is posed in the form 


$ = 


1/3 


[+ 1< 1 * * * * * 7 ’ t) cos tot + x^(r*, t) sin cot J 


+ e £i(j J[r) + 4 * 2 ( 7 ) cos 2u>t + x^fr) sin 2tot| 
+ « [+ 3 <?> cos 3tot + x^(~rj sin 3totj , 


(A. 18) 


1 

In the preliminary efforts to find a solution that was valid through 
third order terms, a solution which included only the Jj modes was 
first posed. This, of course, is sufficiently general to obtain an exact 
solution for the linear problem, but is not adequate for the nonlinear 

problem. This work clearly showed that the fluid problem is related 
to the Duffing problem, in that two second order differential equations 

to be satisfied by the generalized coordinate Ajj and Bn were of the 

Duffing type and might be described as Duffing's equations generalized 

to two dimensions. However, more importantly, this work demonstrated 

that other J m modes are excited through nonlinear coupling of the free- 
surface waves. This coupling involved a coupling between the Jj modes 
and the other J m modes that started with third- order terms (only J Q and 
J^)- It was also noted that the J Q and modes had a steady- state 
response that was second harmonic in time. These results were used 
as a guide in posing the steady state solution of boundary value problem 
[Equations (A. 12) and (A. 13)] in the form given in Equation (A. 18). 
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where 4< n and x n satisfy Laplace's equation and only ijjj and Xj depend 
upon time. The notation!* means dependance upon r, 9, and z. 


In this development, the arbitrary potential functions, 4* n and 
are chosen to identically satisfy Equation (A. 14). The time varying 
generalized coordinate associated with each of the potential functions 
is selected to approximately satisfy free-surface condition Equa- 
tion (A. 15). The approximate satisfaction of Equation (A. 15) is carried 

out by first introducing Equation (A. 18) into (A. 15) and then requiring 

1/3 2/3 

that the coefficients of € , e ' , and e vanish for all time in a 

Rayleigh- Ritz sense. In the following work, it is noted that the coeffi- 

l/3 

cient of « contains sin cut and cos cut terms; that the coefficient of 
2/3 

f contains sin 2cut, cos 2cot terms; and terms independent of cut; and 

finally, that the coefficient of e contains sin cut and cos cut terms and 

also higher harmonic terms. The approximation imposed is that these 

l/3 2/3 

terms all vanish identically for the « and € terms and the first 
harmonic terms for € . 


When these steps are carried out, the transformations 


1 2/3 2 2. , 2/3, /A 1Q , 

T = -^ € CUt>Pn=Cu(l-V€ ) (A. 19) 

are introduced where the lowest natural frequency of small, free 
surface oscillations, p^j, i® given by 


Pll 



tanh ( X jh) 


(A. 20) 


and where corresponds to the first zero of J^' (X^a). Introducing 

(A. 18) into (A. 15), using (A. 19), and then setting the coefficient of the 
1/3 

c terms equal to zero yields 
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(g i|/ lz - Pj ] 2 iJ/j] COS wt + ^gXj z - Pj 1 2 x 1 j sin cot = 0 on z - 0 . A. 21) 

For these first-order terms to vanish for all time, it is 
necessary that 


* 


'll 


1 z 


tj ^ and Xj 


'll 

x, 

g 1 


(A. 22) 


Equation (A. 21) will be satisfied identically by choosing 

cosh (\ (z + hfi 

+1 ° [ f l< T > c ° s 6 + f 3< T) 3in e >l< X ll r » ■ bosK i ^,M ■ ■ 


cosh[\j,(z + h)l 

X 1 = [ f 2 (T) COS 9 +f 4 (T) 8in e J J l (X ll r) cosh (k u h) 


V (A. 23) 


regardless of the values of the generalized coordinates f., providing 


Pjl " ^ll a ll 9 a l 1 ~ tanh(Xjjh) 


(A. 24) 


2/3 


For the second- order terms to vanish, the coefficient of c 
must vanish. The resulting equation contains second harmonic terms 
in wt and constant terms only. Setting the constant term and the 
cos (2u>t) and sin (2cot) terms equal to zero and using the relations 


x l*lz = x lzh and X l*lzz = X lzz*l 


(A. 25) 


obtained from Equations (A. 22), reduces the three equations to 


*^0z 


(A. 26a) 
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2 ( 1 ^ a ll ^ 2 A 

4p 11^2‘ g ^2z = 2p ll\^ x lr^lr + ~2 X 1 © ^ 1 0 + 2 K 1 1 X 1 ^ 1 )• 


(A. 26b) 


4p l 1 x 2 ‘ gx 2z = P 


f 2 2 . 1 2 1 

11 X lr " +lr + “2 X 10 _ T+li 

L r r 

„ 2 (-.* ■ *. 2 )] • 


3a i i - 1 

+ i_L \ 


(A. 26c) 


The arbitrary functions 4 ^, 4 ^, and are selected to satisfy 
Equations (A.26a,b, c). If 4 >q is taken to be a constant, Equation (A. 26a) 
will be satisfied identically. If 4 ^ an ^ are taken to be 


00 


r * ... cosh[\ 0n (z + hj 

^2 “ Z_ A 0n J 0* X 0n r * cosh (\ Q K) 

n=l VI1 


00 


>v coshQx, (z + hjj 

+ L (A 2n C ° S 26 + B 2n Sin 20)J 2 (X 2 n r) cosh (X^T) (A ' 27a) 

n=l 


and 


= Z C 0n J 0 (X 0n r) co?MX' n M 


n=l 


cosh^Xg^z + h)J 

— I 
On 


& A A coshfx., (z + h)l 

+ ) (C 2n co, 26 + D 2n sin Z W { X 

n= 1 ' 2n & * * 9 


where 


(A. 27b) 


J i< X 0„*> = J 2<^n a > = 0 • 


(A. 28) 


then Equations (A. 26b) and (A. 26c) can be satisfied by choosing the 
suitable generalized coordinates in 4 ^ and x^. These generalized 
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coordinates can be expressed in terms of f ^ . f^, iy a nd ^ by intro- 
ducing Equations (A. 23) and (A.27a,b) into (A.26b,c) and applying the 
Fourier-Bessel techniques in which the following orthogonality relations 
are used: 


f rj (\ n r ) J (\ n r )dr 
0 Om 0 On 
Jo 


f IJ 2 (X 2m r)J 2< X 2n r > dr = 

J o 


. m 4 n 


a m = n 


sT J 0 ' On 


f. 


2\ 


2n 


(A. 29a) 


m * n 


x 2 2 A 
X_ a - 4 

— r— J 2 ^ 2n a) - m = n 


(A. 29b) 


In this manner, the generalized coordinates of the Jq and modes 
can be expressed in terms of the E. This process yields 


A n = (f.f, + f,f J , C n 

On On' 1 2 3 4 On 


n 0„7( f 2 2 +, 4 2 - £ l 2 - £ 3 2 ) ■] 


*2„ - n 2„< £ l £ 2 - £ 3 £ 4> ■ £ 2„ = n 2 n l(‘l + l 3 ' l l ' l i) ’ 

Kn - n 2„< £ l £ 4 + £ 2 £ 3> • °2n = n 2n( £ 2 £ 4 ‘ £ 1 £ 3 ) ' 


(A. 30) 
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where 


Q 


T ^ X T ^ i. t/ I ® 

A 01 02 + 03 


On 


(a 2 2\ a T 2 

("Pll -Ponj^T/o 


< x o„*> 




T n T n + K T n 

hi ' hz +k o a 23 




2n 


( 4 Pll 2 - Pz n 2 ) 


x 7 T . 

X 2n a - 4 


r 

2n P 11 


2X-, p 


J 2 < X 2n a > 


J 


(A. 31) 


and 


n 


ql 


n 


q2 


T n 

I Q = 

q3 


K 


f Uv-il 


2 

du , q = 0, 2 , 


r X 1 I a 


r X 1 I a 


i J q(x77 U ) J l 2 < u ) du * q - 0, 2 , 


11 


uJj t »^ — u 1 Jj^(u)du , q = 0, 2 , 


\ (A. 32) 


q V 11 


3a ll " 1 


The coefficient of € yields the third-order terms. Third-order 
terms arise from each of the boundary value terms, Bj, B^, and B^. 
For the third-order terms, it is only required that the first harmonic 
terms vanish. The first harmonic terms from B^ are 
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’ll 


f dx l 
l dT 


- _ r cos 0j cos wt 


(“3T +vx l) 8in ut ] 


(A. 33a) 


and the first harmonic terms from are 


’ll 


' 1 'lr x 2r- x ]r ,t '2r + -7 ( ' t, ie x 2e- x 19' t, 2e ) - X H (“ll - *) <*1+2 ‘ Vz> 

r 

- X 1 l a l l'^Zz ‘ x 1^2z> ^l^zz ' x l+2zz>] cos “* ' Pi l^lr^r + + 1 r+2r 
+ “Z ,x i9 x 29 + ' |l i9'W + K i i(“i 1 ' 1 )< x i 1 ‘ 2 + ■h+2> ' X u a n (x i x 2z 


+ +1^ 2 z + ? ( x i x 22 z+ + 


l't'Zzz 1 ] 31 


in wt . 


(A. 33b) 


Similarly, there are first harmonic terms from which are 
lengthly expressions and are not presented because it is the integrals 
of Bj, B^, and B^ that are required in the Rayleigh-Ritz procedure. 
Introducing Equations (A. 23) into the expressions for Bj, B^, and B^ 
and using Equations (A. 30) to evaluate the contribution from B^ results 
in the following integrals: 

-a f 2 tt 


n CTT 

Bj cos 6 J^XjjrJr dr d6 

3 


= Trp 


11 


^ 1 1 a ” * 2 (^2 ^ 
T~ J 1 (X ll a) \dT ‘ vf lJ 


2 \ 


11 


T J l (x il a * 


L ll 


cos wt 


- Trp 


11 


. 2 2 
K 11 a 


2X. 


11 


“ J i 2(x ii a) (-^r + vf 2) 


sin wt , 


(A. 34a) 
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a fZ-n 



B 


o J o 


j sin 0 Jj(X^r) r dr d0 


^11 


; 2 2 . 

X n a - i 2 


2X 


11 


f d{ 4 ' 

2 J i (x n a) (^r ‘ Vf 3y 


cos u)t (A. 34b) 


^11 


■a :'2 tt 


X, , 2 a 2 - 1 2 /«, 

L 2 *il 


sin wt , 


I I cos 0 J^X^rJr dr d0 

J o J o 

= " Pn [ f i( f i 
+ f 4 (f 2 f 3 - fl £ 4 )G 2 J 


2 + f . 2 + f . 2 + f/JSj 


2 ’ A 3 ’ 4 
cos cat 


+ "i>n[ f 2 ( f i 2 * l z + J 3 2 + */ 

~ f 3 (f 2 f 3 " f i f 4 )G 2 ]sinut , 


A 


(A. 34c) 


a rZ.tr 



o J o 


B^ sin 0 J^(X^r)r dr d0 


[ £ 3 ( f . 2 


- * 3 ^i + f 2 + f 3 + f 4 


ft 


" f 2^ f 2 f 3 " f 1 f 4^ G zJ COS Wt 
+ "I’ll Q 4 ( f l 2+ l Z + { 3 + f /)°l 

+ f l^ f 2 f 3 “ f l f 4 ^ G 2 ^| Sin Wt * 


(A. 34d) 
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where 


a /tZtt 



o ' o 


cos 0 dr d0 


7 p n 2 f i( f i 


+ W 3 - f 


2 + f/ + f 3 2 + f/jHj 

cos cot 


1 f 4 ) H 2 J c °s cot -^Pn 2 ^! 

3 


2 +f 2 +f 2 
2 3 


+ f / A ! - W3- t . f 4»«2| Sin “ t ’ 


a /'2'rr 



o so 


sin 6 Jj(\^r)r dr d0 

' ■J p n 2 [ f 3( f i 
- £ 2 (£ 2 £ 3 • r i f 4 ) “I] cos 
-7 p ll 2 [ f 4( f l 2 + J 2 2 + f 3 2t£ 4 2 )»l 
+ £ 1< £ 2 £ 3 - £ 1 £ 4>« 2 ] si 


2 + f 2 2 + I 3 2 4 f/)A, 


sin cot 


(A. 34e) 


(A. 34£) 
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G, 
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n=l v l - 


„ . . 1.2 

Q 0n 1 1 On 1 1 2 On 


+ \ 


11 


7 n 2 n J 24 


+ X. 


11 


(* ’ a n 2 )] 

11 + z[ a ll a 

(> - “n 2 )] °Znh 3 j • 


n 0n I 03 " ^0n I 04 " n 2n*22 


XX 1 \ 2 

2n*ll*2n " 1 K Zn 


(A- 35a) 
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(A. 35b) 


(A. 36) 
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(A. 36) 
(Cont. ) 
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Two ordinary differential equations are obtained by introducing 
the potential function Equations (A. 18) into (A. 15) and then carrying 
out the Rayleigh- Ritz procedure by first multiplying by Jj(Xjjr) 
cos 0 rdr d0 and integrating over 0 -s r a and 0 < 9 < 2tr, then 
multiplying by J^(\^r) sin 0r dr d0 and integrating throughout the 
same region. Then, the first harmonic terms in w are set equal to 
zero in both of these equations, producing the following four first- 
order nonlinear differential equations: 


= 

dT 

" vf 2 

- K i f a( f i 2 + ‘z + h 2 + 

) + K 2 f 3^ f 2 f 3 

- f l f 4> • 1 

d f 2 

F 1 + 

v£ l + K l £ l( £ l + £ 2 + £ 3 

+ f 4 ) + K 2 f 4^ f 2 f 3 ' f l f 4* ’ 

f3 = 

dT 

- vf 4 

- K l f 4 ( f l 2 + f 2 2 + f 3 2 + £ - 

4 ) " K 2 f l* f 2 £ 3 

' f l f 4> ’ 

^4 = 
dT 

vf^ + 

K l £ 3( £ l 2+£ 2 + £ 3 + £ 4 ) 

^ " K 2 f 2^ f 2 f 3 ‘ 

f l f 4> ' J 

where 






K 1 = K 10 + AK 1 


K 2 = K 20 + AK 2 


F, = 


2a 


1 


(A. 37) 


(A. 38a) 
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and 
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-181, + 31, + 61 . + 9I C - 12L + 6 l_ 
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L 11 
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l ll 
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The form of Equations (A. 37) is identical with that of the set of 

5 

differential equations derived by Miles for the undamped spherical 
pendulum. Therefore, the steady-state solutions of Equations (A. 37) 
are the same as the two obtained by Miles for the spherical pendulum. 

These two steady- state harmonic solutions correspond to the zeros 
of dT/dT for i = 1, 2, 3, 4. The first solution, called planar motion, 
is a steady- state fluid motion with a constant peak wave height and a 
single, stationary nodal dimeter perpendicular to the direction of ex- 
citation. The second solution, called nonplanar motion, is a steady- 
state fluid motion with a constant peak wave height and a single nodal 
diameter that rotates at a constant rate around the container. For the 
planar motion, the solution is 


= 7 • f 2 = f 3 = f 4 = 0 ’ 


(A. 39a) 



40 


with 


„ - 1 T 2 
v = - f jY - KjY . 


(A. 39b) 


For the nonplanar motion, the solution is 


f 


1 




z , F i 

V +K - 2 


-1 

Y 


(A. 40a) 


with 


v = - K 3 y _1 + K 4 y 2 , (A. 40b) 

where 

■ TJ F 1 • K 4 = K 2 - 2K 1 • (A ' 41) 

3 

This solution is real, and hence exists, for \ =» 0 when Y + 0 

3 i c. 

and for y< 0 when y + 0. 

4. Stability of Steady-State Harmonic Solutions 

Stability of the steady- state harmonic solutions obtained in the 
preceding section will now be investigated. A particular steady- state 
harmonic solution of Equations (A. 37) is denoted by the superscript (o), 
and the stability of such solutions is investigated by imposing a small 
perturbation from this steady- state solution and examining the subsequent 
motion. If the motion following the perturbation decreases with time, 
the solution is called stable; if the motion increases with time, the 
solution is called unstable. Mathematically, set 


^(T) = f. (0) + 


\r 


c. e 
i 


(A. 42) 
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where £.^°^ are constants corresponding to the steady- state amplitudes 
of the harmonic solutions of Equations (A. 37) and the c. are assumed to 

i \ 1 

be small. Stable E' solutions correspond to values of X with negative 
real part, and unstable solutions correspond to correspond to values 
of X with positive real part. Introducing Equation (A. 42) into Equations 
(A. 37) neglecting products of c^, and imposing the condition that the 
f^°^ are solutions of Equations (A. 37) leads to the following set of 
homogeneous algebraic equations: 


d ll + x 

d 12 

d l 3 

d 14 


r '' 

c i 


r <f 

d 21 

d 22" X 

d 23 

d 24 

< 

C 2 

>-< 

0 

d 31 

d 32 

d 33 + \ 

d 34 


C 3 


0 

_ d 41 

d 42 

d 43 

d 44 • \ 




0 

V J 


(A. 43) 
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where 


2K 1 f 1 ( °^° ) + K 2 f 3 (o) f 4 (o) , 


v + K.f/ 0 ' 2 ♦ 3f 2 <°> 2 ♦ f,«» 2 + f 4 «°» 2 ]- K 2 f 3 <°> 2 , 

+ K 2 jfj ( °>/ 4 < o) . Zf 2 (o) f3 <o) ] , 
2K 1 f 2 (o) f^ o) + K 2 f 1 ( °^° ) , 


v + K 




l ll ' 


= 2K 1 f 1 (o) f 3 (o) + K 2 f 2 (o) f 4 (o) , 


2K 1 f 1 (o) f 4 ( ° ) + K 2 |f^ 0 ^f 3 (0) - uM'to] , 


43 


d 31 " d 24 ' 


‘32 


*33 


‘34 


= d 


14 ’ 


2K l'3 0)£ 4 0) + K 2 f l <0)f 2 0) • 

r f (o. 2 + f <o . 2 1 f «o»% 3f (o» 2 ] . K 2 r <°« 2 . 


= V + K 


d 41 “ d 23 ’ 


‘42 


‘43 


= d 


13 ' 


* 


v + K,|£.<°> 2 + £ 2 < 0 > 2 + 3f 3 < 0 > 2 


+ £ 4 «°> 2 ]- 


K 2 f 2 


(o) 


d 44 " d 33 ‘ 


V (A. 44) 
' (Cont. ) 


J 
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Since Equation (A. 43) is a homogeneous set of equations in c., 
nontrivial solutions will exist only if the determinant of the coefficients 
in zero. Setting this determinant equal to zero yields the allowable 
values of X. The question of stability reduces to an examination of the 
roots of the resulting equation in X. 

a. Planar Motion Solution 


The steady- state planar motion solution found in this section 
was given as 


f 


(o) 

1 


= X 




(A. 45a) 


where 


v = -FjY - 1 - K lV 2 . (A. 45b) 

Introducing Equations (A. 45a) into Equation (A. 43), setting the 
determinant of the coefficients of c. equal to zero, and expanding 
results in 


where 


X 4 + (Mj + M^JX 2 + Mj - 0 , 


(A. 46) 


Mi = 


M. 


_ 2 



v + KjY 

0 


0 

v + 3 K } y 2 


v + Kj Y 2 

0 


0 

v + (Kj - K 2 )y 2 

J 


> (A. 47) 
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The roots of Equation (A. 46) are 


= -M, - - (v + K lY 2 )(v + 3 KjY 2 ) - -F 1 y" 1 (f,y" 1 - 2K lY 2 )< A 48a ) 


and 


V ' - M 2 


= -(‘' + K l' 2 ) 


v + (Kj - K 2 ) Y 


= -Fiy' 1 (f iY _1 + k zY 2 ) . (A .48b) 


The boundary between stable and unstable planar motion cor- 
responds to = 0 and = 0. From = 0, it is found that: 


vl/3 


Y = ±oo and Y 


2K, 


(A.49a) 


Later, it is seen that Kj =*• 0, so that from Equation (A. 45b) these 
values correspond to 


v = -oo and v 


iZ/3 


-3K 


1 


1\ 2K 


1 , 


F rom X- 2 


2 


= 0 , 


y = ±co and y 



(A. 49b) 


(A. 50a) 


are obtained. 

When Equation (A. 45b) is used with Kj =*■ 0, it is found that 
these values correspond to 



2/3 


v 


-oo and v = (K 2 - K x ) 


(A. 50b) 
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Figure 2 summarizes the regions of stable and unstable planar 
motion. 

b. Nonplanar Motion Solution 

The steady- state nonplanar motion found in this section was 

given as 


f 


(°) 

1 


Y- 


f (°) _ f (°) 
2 " 3 




2 

Y 



(A. 51a) 


where 

v = - K 3 y - 1 + K 4 y 2 (A. 51b) 

2 

The requirement that t, be real demands that f, =■* 0, which, in turn, 
indicates that this solution is not valid when y lies in the range 


-Y 3 0 • (A. 52) 

^2 


Introducing Equations (A. 51a) into Equation (A. 43), setting 
the determinant of the coefficients of c^ equal to zero, and expanding 
results in 


\ 4 + 


(m 3 + 


m 4 )\ 2 + 


m 5 m 6 = 


(A. 53) 
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where 



(A. 54) 


and 

d l 2 = v + Kjfv 2 + i Z ) . 

d 13 = K 2* l. 

d 2i = v + 3KjY 2 + (Kj - K z )i Z , 

d 24 = 2 < K l " K 2>^ • 

d 34 = v + 3K i^ 2 + < K i - K 2 )Y 2 . 


Introducing Equation (A. 51b) into Equations (A. 55) and ex- 
panding Equations (A. 54) yields 





REFERENCES 


Eulitz, W.R. and Glaser, R.F., "Comparative Experimental and 
Theoretical Considerations on the Mechanisms of Fluid 
Oscillations in Cylindrical Containers." Report MTP-M-S and 
M-61, Army Ballistic Missile Agency, Huntsville, Ala., 

29 May 1961. 

Berlot, R.R., "Production of Rotation in a Confined Liquid Through 
Translation Motion of the Boundaries," J. Appl. Mech., 
26:513-516, Dec. 1959. 

Milne-Thomson, T.M. Theoretical Hydrodynamics , 3rd ed., 
Macmillan, London, 1950. ~ 

British Association for the Advancement of Science, Bessel 
Functions , University Press, Cambridge, Eng., 1950-58, 

2 vols. (its Mathematical Tables , vol. 6 and 10) Part I, 
Functions of Orders Zero and Unity; Part II, Functions of 
Positive Integer Order. 

Miles, J. W., "Stability of Forced Oscillations of A Spherical 
Pendulum," Quart. Appl. Math, 20:21-32, April 1962. 














Table 2. Computed Values of Parameters 


50 



r- 

iTi 

m 

o 

vO 



n- 

00 

o 

'O 

eg 
n i 

CO 


ON 

00 

in 

O 

sO 

o 


vO 

" r\j 

IT) 

m 

m 

o 

in 


00 


t— ( 


o 


O 

o 

o 

o 

o 


o 

© 

o 

o 

d 






eg 


sO 

t- 


00 

eg 


eg 

in 

eg 


in 


r- 

eg 

rj 

m 

00 

fH 

'O 

sO 

o 


CO 

d rj 


h- 

I s - 

o' 

vO 

>-H 

r\J 

co 

H 

o 

O 


O 

O 

o 

o 

O 


o 

d 

o 

o 

o 






00 




eg 

^H 

•—< 




co 

CO 

o 



M 

0 s 

CO 

O' 



I s - 


n- 

o 



vO 

r- 

ro 

o 

c 5 

o 

'T 

m 


o 

l-H 



o 

O 

o 




o 

o 

o 



o 

o 

d 

o 





o 

<-H 

r- 

■ 



o 


-H 



eg 

,H 

l-H 


B 


CO 

00 

r- 

o 


^_i 

O 

eg 


CO 


00 

m 

00 

in 

o 

BSiB 

o 

O' 

o 

o 

o 


pg 

vO 

o 

o 

o 


V 

o 

i 

o 

o 

o 






CO 




o 


CO 



eg 


m 

0 s 


CO 

C0 

eg 

i-H 

CO 


00 

eg 


eg 

m 

eg 

r- 

NO 

m 

O' 

00 

d o 

o 

00 


1-4 

o 

►H 

00 

eg 

o 

o 

o 


rg 

o 

o 

o 

o 


o 

o 

d 

i 

o 

o 

■ 




o 





o 

0 s 

o 




00 

«H 

eg 

: ' B 


CO 

r- 

o 

O 

'' B 

(M 

r- 

vO 


t- 


o 

o 

c- 

CO 

-H 


00 

eg 

CO 

O 

o 


CO 

o 

o 

O 

o 


rg 


o 

o 

o 


—* 

o 

o 

o 

o 


o 

d 

o 

o 

o 




00 

r- 

o 

00 

eg 

— < 


m 

o 

CO 

r- 

00 


O' 

o 

CO 

O' 

00 

O 

o 

m 

o 

eg 

m 

CO 

rg 

eg 

o 

o 

r- 

in 

CO 


i-H 

—< 


o 

yp 

^*4 

eg 

CO 

o 

r- 

o 

*-» 

o 

o 

o 

o 

O 


o 

d 

i 

o 

o 

o 

o 

o 

o 

o 



eg 

CO 


in 

'•O 

r- 

00 

0 s 



in 



m 



i 

n 


i 



o 

t 


o 



—H 

o 


i-H 









X 



X 


d 

00 

X 





r- 

vO 

CO 

o 



eg 

SO 

m 




m 

o 

CO 

o 



00 

r- 

eg 

o 




co 

o 




o 


CO* 

o 



in 

in 





i 

o 

o 





<— i 





o 

X 

X 




d 

CO 

CO 





00 

00 






00 





O' 

f- 





CO 

r* 





o 

— 









r- 





sO 

eg 




cH 

O' 

CO 


00 

m 

o 

o 

«-H 


CO 

n- 

o 




NO 

O' 

g< 

eg 


eg 

O' 

i-H 

o 


o 


in 

O' 

i-H 

O 

o 

C5 

—H 

o 

o 

o 

o 


o 

o 

o 

o 

o 


o 

o 

o 

o 

o 


d 

0 

1 

o 

d 

o 

i 









CO 


co 



o 

vO 

r- 

r- 

CO 


m 

in 

eg 

O' 

sO 


o 

eg 

r- 

CO 

r- 

d 

g* 

-H 

o 

CO 

eg 

o 

o 

eg 

in 


o 

a 

H 

o 

t-4 

o 

o 


o 

o 

o 

o 

o 


o 

o 

o 

o 

o 


• 

• 

• 

f 

ft 


o 

o 

o 

1 

o 

1 

o 


00 

vO 

00 

vO 

o 


m 

m 

t- 

o 

O' 


o 

i-H 

yp 

00 

H 


O' 

4* 

eg 


r- 

. g* 

o 

O 

eg 

o 

eg 

d rg 

in 

00 

CO 

eg 

i-H 

►H 

o 

o 

o 

o 

o 


d 

0 

1 

o 

o 

o 

i 



O' 

eg 

eg 

00 



r- 

e- 

'O 





i-H 

sO 



* 

yp 

vO 

r- 

l-H 

— CO 

i ** 

00 

i-H 

sO 

eg 

0 rg 

o 

00 

i-H 

m 

r- 

F-H 

n- 

o 

o 

o 

o 


o 

o 

o 

o 

o 


• 

ft 

• 

• 

• 



o 

o 

l 

o 

1 

o 

1 





• 


d 

- 

eg 

CO 

B 

m 


i 


Note: K = 0.97624 . K = 0.96209 x 10 























51 


Table 3. Domain of Stable and Unstable 
Nonplanar Harmonic Motions 


-oo «< v ■< -0. 03027 
0 <y «= 64. 47 
117.6-= U|-=oo 


unstable, since Re(X) is positive 


-0. 03027 -=v<oo 
64. 47 -= y <oo 
117.6 1 1, | -= oo 

0. 06459 ■*= v -=oo 
-oo «=• y -c -85. 41 
0 «| ^1 <00 

oo <v<0. 06459 
-85. 41 -= y <0 


stable, 


since Re(X) is negative 


unstable, since Re(\) is positive 


} 


solution does not exist, since £ «= 0 
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Figure 1. Coordinate Systems 
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STABLE MOTION 
UNSTABLE MOTION 


55 



X 


tion Amplitud 






56 



Figure 5. Nonplanar Motion Amplitude/Frequency Relation, 
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Figure 7. Unstable Regions for Planar and Nonplanar 
Harmonic Motions 




igure 8. Planar Motion Instability Region as a Function of 
Nondimensional Driving Frequency and Amplitude 
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Figure 9. Slosh Test Facility With One-Degree-of-F reedom 
Tank Platform Laterally Oscillated with Scoth- 
Yoke Sinusoidal-Drive Mechanism 
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Figure 13. Comparison of Theoretical and Experimental 
Planar and Nonplanar Instability Regions 


NASA-Langley, 1963 D— X8T0 




